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Abstract 

Intrinsic equations of motion of superconducting cosmic string may 
admit solutions in the shock-wave form that implies discontinuity of 
the current term \- The hypersurface of discontinuity propagates at 
finite velocity determined by finite increment A% = \+ ~ X-- The 
current increases x+ > X- m stable shocks but transition between 
spacelike (x > 0) and timelike (x < 0) currents is impossible. 

1 Introduction 

Cosmic strings are 2-dimensional topological defects that may be formed at a 
phase transition in the early universe [TJ [2]. Although their direct detection 
is still impossible with modern technical equipment, they are believed to be 
responsible for several astrophysical phenomena associated with gravitational 
lensing, gravitational waves, particle acceleration and gamma-ray bursts [3j 
Sj, and their theoretical study can give better understanding of the global 
structure of the universe. 

The dynamics cosmic strings is determined by the Lagrangian A, also 
called as equation of state (EOS). The Lagrangian of the Goto-Nambu 
model A = —m 2 depends on the only mass parameter m. Cosmic strings 
may have non-trivial internal structure, giving rise to additional degrees of 
freedom associated, for example, with transversal oscillations of the world- 
sheet, or actual particles (bosons or fermions) trapped in the defect core 



[5] . Such current-carrying or " superconducting" strings are described by the 
Lagrangian A = A(m,x) which depends on the magnitude of the current 
X = —h ab d<f)/da a d<j)/dcr b that includes the surface metric h ab with respect 
to local coordinates o a within the string worldsheet, and the scalar field 
coupled to the string. 

There are only few explicit analytical models of superconducting strings. 
The linear EOS [6] A = — m? — x/2 is applied to the strings that carry 
fermionic currents. Bosonic current carriers are described by much more 
complicated EOS, namely, the "transonic" model [TJ, [8J, [9] 

A = —my/ m 2 + x (1) 

the "polynomial" model [TO] : 

A = - m2 -f( 1 + ,|) < 2 > 

the "rational" model pH [12] : 

A -- m2 -l( 1 + ^f)" < 3 > 

and the " logarithmic" model [TOl [13] : 

A — ! -y'"( 1 + 4) (4) 

where additional mass parameter is introduced. 

In contrast to the Goto-Nambu model, the superconducting string models 
admit a large variety of stable loop solutions (or, vortons) and their inter- 
commutation may frequently take place [HI ITS] . The problem of vorton 
stability and evolution of superconducting string networks attract the inter- 
est of researchers because such string loops can accumulate significant part 
of the universe mass and the expected dominant effect of their reconnection 
is particle radiation (when some particles are expelled away). 

The equations of motion of superconducting strings admit solutions in the 
form of infinitesimal perturbations of two types [16]: extrinsic (transverse) 
perturbations of the worldsheet that are oscillations of the string geometry, 
and intrinsic (longitudinal) perturbations within the worldsheet that are os- 
cillations of the current x- However, this equilibrium, or so-called "elastic" 
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description is admitted until small transverse oscillations grow up and form 
abrupt changes in the geometry (kinks), or when small longitudinal oscil- 
lations are accumulated into finite-amplitude discontinuities of the current 
(shock waves) [T71[T8]. Such strings loops may tend to fold on themselves and 
make contact points of self-intersection so that it will be energetically favored 
to some of the trapped particles to move out of the string, thus, generating 
emission, associated with possible visible consequences [U El El S] • 

The dynamical evolution of a fermionic current carrier can develop only 
kinks, while a bosonic current carrier can develop kinks at time-like (or "elec- 
tric") currents x < and shocks at space-like (or "magnetic") currents x > 
[TTl [T8] . In the latter case nonlinear effects become dominant rapidly in lon- 
gitudinal perturbations and a stable discontinuity appears: the magnitude 
X becomes discontinuous and the total energy is no longer conserved. This 
singular behavior of superconducting strings cannot be assigned to some pro- 
jection effect because it is seen in simple analytical analysis [19J as well as it 
is observed in numerical simulations performed both in two dimensions [T7| 
and in three dimensions [18] . An evident hint to a shock wave can be also 
recognized in numerical simulation of a 2-dimensional vorton in the magnetic 
regime (Figure 7 in Ref. [20]). 

In the present paper we complete the analytical analysis of shocks in 
superconducting cosmic strings. Particular interest is focused on the behavior 
of the bosonic current, the shock velocity, the transition between electric and 
magnetic regimes. 

2 Small perturbations and shocks 

The physical behavior of superconducting strings is described by the equa- 
tions of motion. Their projection along the string worldsheet yields the 
'intrinsic' equations or the conservation laws j!6] 



for = /iv^ and j M = nv M where tensor rf„ = v v v il — u v u^ is composed of 
mutually orthogonal unit vectors (u^u^ = —1 = —v^v^ and u^v^ = 0), and 
dynamical parameters [UJ 



(5) 




n 2 = K 2 x 
n 2 = -x 



X>0 
X<0 



(6) 
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are determined by function [TD] 

The longitudinal or 'sound' perturbations run within the worldsheet at 
velocity [E] c 2 = c\ = (n/fi) dfi/dn that in the light of @-@ is written in 
the form 

9 ( Y dK\ Tl 

c 2 = 1 + 2-^— (8) 



Kd X 

where the upper (lower) signs correspond to x > (x < 0). It is always 
K < 1 at x < 0, K (0) = 1 and K > 1 at x > because dK/dx > for all 
known models (fl|-(j3j • 

According to the linear EOS of fermionic cosmic strings A = —m 2 — x/2, 
the sound speed ^ coincides with the speed of light c = cu g ht = 1. The 
bosonic string models ([l])-(j3]) are endowed with finite sound speed dispersion 
dc/dx 7^ that provides a possibility for development of shock waves [19] . 
As a matter of fact, the fermionic current-carrying strings can develop only 
kinks, while the bosonic strings admit both kinks and shocks [T7J [IB"] . 

The front of a shock wave is a hypersurface in 4-dimensional space whose 
equation is given in the form <p (x) = 0. It implies that + = 0_ for any two 
points in the hypersurface but, in general, <p + — </>_ ^ for arbitrary two 
points beyond the hypersureface, particularly, when "— " and "+" label the 
states before and behind the front. A stable shock front is associated with 
characteristic space-like vector 

which is the same before and behind the front X +u = = X u . It can be 
presented in the basis (u u ,v v ) in the following form 

_W-U„ V + V- V _W + U +U + V +I/ 
v - / -, 2 - 7j 2 l iU J 

v 1 - w - v 1 - ^+ 

where parameters w_ and w + play the role of velocities whose essence is 
clarified below. 

All dynamical variables in the intrinsic equations of motion ^ are de- 
pendent on function <p(x), and, crossing the discontinuity, we have 



v+ l y^ = v-^uf-+[ v ;v v r] (ii) 
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where square brackets imply a jump of parameters across the discontinuity. 
In the light of (|5]), equation (11) implies 



% 



Df 



d(f) dx v 



and at finite + — 0_ 
that implies X +fl ^i + vt 



f- 



7^ it yields relation at the front A+^j' 



(12) 



A_ M/ u_f^ and A +/J n + w+ 



A_„n_M^. Taking (10) 



we find the shock parameters 



w 



n\ n\ - fj?_ 
fi\ n 2 + — n 2 _ 



an 



n 2 _ /i 2 . 



M 2 _ 



jjL 2 _ n\ 



n 



(13) 



Consider a discontinuity in the the laboratory reference frame (see scheme 
in Fig. [IJi) and, altogether, in the reference frame, co-moving its front (see 
Fig. [TJd). The front is at rest and the flow before it has finite velocity, say, 
W-, while the flow behind the front has velocity w+. When we turn to the 
laboratory reference frame (Fig. [IJi), the shock- wave propagates along the 
string at velocity D- = —W-, while the velocity behind the front D + = 
(w + — to_)/(l — w_w + ) is finite if w + ^ w. 
w + — > c in the sound limit n + — > n_ and fi + — > 
and D + = corresponds to a sound wave ^ because no material flow is 
generated behind the shock. 

Taking /z and n from equation (|6]), we determine the shock wave velocity 
in the magnetic regime (x > 0): 



Formula (13) yields — > 
and /i + — > Hence, D_ = —c 



w_ 



K 2 



Ka 



x+ - X- 



■Kix + -Ki X - 



(14) 



and in the electric regime (% < 0): 

K 2 



ur 



1 K 2 +X+ -K 2 _x- 



x+ - X- 



(15) 



3 Growth of current 

Existence of a stable shock is possible when the evolutionary condition is 
satisfied [22J E3]. For the shock waves in strings it implies [T9] : 



W- > c_ 



w + < c + 



(16) 
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that coincides with the same criterion for the ordinary shock waves. However, 
this fact can be predicted immediately by means of very easy analysis. The 
number of intrinsic equations of motion ^ for = fxv M and j M = nv^ is two, 
while there is only one material parameter x~- hence, the shock in a string 
has exactly one degree of freedom. The shock waves in continuous medium 
are described by three independent equations of motion [221 while there 
are two material parameters (the pressure and entropy), and, again, there is 
only one degree of freedom. It implies that the shocks in fluids and cosmic 



strings obey the same constraint (16), imposed on their velocities. 

Now consider small-amplitude shock waves. If increment Ax = X~\ 
is small with respect to |x±|, formulas (14)- (15) can be expanded in a series 



X- 



of Ax- According to (14)- (15), we have 



w 



1 - 



Q 



K 



X- 



> 



(17) 



w_ 



1 



c 



, §2X-^X 



x-<o 



where 



q_ 3 ( dK Y -K— 

\dx J dx 2 



(18) 



(19) 



and K is determined by (ffl). Quantity Q is positive for all known analytic 
EOSs, namely, equation (|TJ yileds Q = mr 2 / (m 2 + x), equation ^ yileds 
Q = 4m^/ (ml — 2%) 4 , equation (JiJ yileds Q = 10 (m 2 + xf / m *> equation 
([I]) yields Q = 3/m^. Hence, in the light of (16), the current increases behind 
the shock wave 

Ax = x+ ~ X- > (20) 

In other words, any stable discontinuity of the current is associated with 
a shock wave propagating in that direction which corresponds to positive 
increment of X- 



Equations (|14j)-(|15j) 
the velocities: 

> 



(16) determine the following relationship between 
in the magnetic regime and c + > 



tw_| > c_ > c + > \w + \ 
W-\ > c_ in the electric regime. Detailed calculation according to 

All EOSs of bosonic strings 



\ w +\ 

formulas (14)-(15) is shown in Fig. |2| and |3 
([!])- ([3]) reveal similar behavior 
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4 Transition between two regimes 



Consider a transition from the electric to magnetic regime. Let us imagine 
it in the form of a joint shock composed of two consequent shocks (Fig. 
Let the first electric shock propagates along the string X- < Xe < and 



velocities W- and w e its determined according to formula (15). Let the second 



magnetic shock propagates immediately after the first shock, x+ > Xm > 



and velocities w + and w m are determined according to formula (14). Both 
shocks merge into a composite shock in the limit \ e — > \ m — > 0. Let the flows 
before and behind the shock are D_ < 1 and D + < D_ < 1, respectively. 
Altogether we determine velocities io_ = —D^ and 



in the co-moving reference frame (Fig. lib). Equations (14)-(15) yield 
— 1 and w m — y — 1, while the flow behind the first electric shock is 

W- — w e 



We. ->■ 



D e = -> 1 (22) 

1 — W P W- 



in the laboratory reference frame (Fig. [Tp.) . The flow behind the second 
magnetic shock is 

D + = W+ ~ Wm -+ 1 (23) 
1 - w m w + 

in the reference frame comoving D e and altogether it is 

D + + D e 

D + = + 6 -> 1 24 
+ 1 + D + D e 1 J 

in the the laboratory reference frame. Therefore, the velocity behind the sec- 
ond shock D + exceeds the shock wave velocity -D_ that implies impossibility 
of a joint shock with a transition from the electric to magnetic regime. 

The limiting chiral regime x = plays the role of a barrier that cannot 
be "tunneled" by a single shock, and the latter cannot be a switcher between 
magnetic and electric regimes of the current. The similar repulsive charac- 
ter of chiral current \ = is observed in the equilibrium solution of the 
field-theoretic vorton model [20] where no transition between the regimes is 
possible: the current vanishes \ — >■ when the loop radius tends to infinity, 
however, the sign of x is n °t changed. The shocks obey the same strange 
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property: it is not able to pass through the chiral point, maintaining con- 
tinuous link between the two regimes. It should be also noted that exact 
chiral limit \ = is no more than a theoretical possibility: as soon as it is 
achieved, the shock becomes asymptotically unstable at finite time interval 
because formulas (14)-(15) give w — c — 1 rather than strict inequality (16) 
which warrants stability at arbitrary time interval. 



5 Conclusion 

The bosonic superconducting strings admit stable discontinuities (shock 
waves) of the current. The shock velocity at space-like currents x > ("mag- 
netic" regime) and time-like currents x < ("electric" regime) is determined 



by formulas (14) and (15), respectively. Small-amplitude shocks are approx- 



imated by formulas (17) and (18). Shock waves of arbitrary magnitude are 
calculated in Fig. [2] and Fig. |3j No shock can provide a transition between 
the electric and magnetic regimes, although electric shock is not forbidden by 



constraint (20). Numerical simulations [TTl IT8] does not reveal existence of 
electric shocks at all, and the physical reason is hidden in their vulnerability 
to transversal oscillations of the string worldsheet [23] that is a problem of 
special investigation beyond the present study. 

The shock propagation implies increase of the current magnitude x+ > X- 



(20). It is similar to increase of the entropy in ordinary shock waves because 
no stable discontinuity is admitted if the current decreases. Condition (}20| 
is equivalent to K + > fT_ for all known models Q-Q. According to (|6j) it 
implies n + > n_ in the magnetic regime. The latter inequality is immediately 
recognized in numerical simulations where a shock structure is associated 
with abrupt jump from lower n_ to higher n + (see especially Fig. 10 in Ref. 

El)- 

In practice, shock waves can appear during non-linear evolution of small 
longitudinal perturbations in the intrinsic equations of motion ([5]). Intercom- 
mutation of two string loops with different currents X- > and x+ > may 
also trigger a shock. Reconnection of two loops with x+ < and X- < 
cannot trigger a shock [TFJ [18]. Reconnection of two loops with opposite 
signs of the current (say, x~ < an d X+ > 0) wm result in appearance of the 
chiral state x = which must remain the same without regard of the radius 
[20]. Although the loop radius must increase and tend to infinity as soon as 
X — > 0, neither smooth nor shock transition from the electric to magnetic 
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regime is possible. We may expect an irreducible point of inflection at x — 
and zero curvature (infinite radius), dividing the electric and magnetic sides 
of the string, while the shock itself can exist at the magnetic side. This prob- 
lem deserves more attention in application to scenarios of possible observable 
phenomena. 

We have considered an ideal stationary problem with no external force. 
The real cosmic string networks may include time-dependent parameters and 
numerical analysis. The analytic: Whenever a shock wave is formed, the 



current is subject to increase (20), and this general law may give help in the 
further theoretical and experimental research of superconducting strings. 

We wish to thank Alexei Starobinsky and Konstantin Stepanyantz for 
support. 
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Figure 1: (a): The shock wave [shaded] propagates at velocity £)_, the flow 
behind it acquires velocity D + . (b): The shock velocities w + and W- are 
determined in the co- moving reference frame, (c): The joint shock transition 
is linking currents X- < an d X+ > 0- 
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Figure 2: The ratio of shock to sound velocity M_ = u> 2 /c 2 (left) and 
M + = w\/c\ (right) vs increment Ax in the magnetic regime of bosonic 
currents at m — m* = 1 at initial x_ = 0.1 are plotted for four models of 
EOS: dot - Q, dash - g, dashdot - Q, solid - Q. 
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Figure 3: The same plots as in Fig. (J2J) but in the electric regime of bosonic 
currents and at initial x_=-0.05. 
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